Phonon dispersions and vibrational properties of monolayer, bilayer, and trilayer 

graphene 
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The phonon dispersions of monolayer and few-layer graphene (AB bilayer, ABA and ABC trilay- 
ers) are investigated using the density-functional perturbation theory (DFPT). Compared with the 
monolayer, the optical phonon E2g mode at F splits into two and three doubly degenerate branches 
for bilayer and trilayer graphene, respectively, due to the weak interlayer coupling. These modes are 
of various symmetry and exhibit different sensitivity to either Raman or infrared (IR) measurements 
(or both). The splitting is found to be 5 cm~^ for bilayer and 2 to 5 cm~^ for trilayer graphene. 
The interlayer coupling is estimated to be about 2 cm~^. We found that the highest optical modes 
at K move up by about 12 cm~^ for bilayer and 18 cm~^ for trilayer relative to monolayer graphene. 
The atomic displacements of these optical eigenmodes are analyzed. 

PACS numbers: 



I. INTRODUCTION 

In recent years, monolayer and few-layer graphene 
have attracted great attention due to the unique proper- 
ties observed experimentally.—'^'^ Many intriguing trans- 
port phenomena, such as balhstic transport at room 
temperature,^'- the anomalous quantum Hall eSect^^ 
and novel many-body couplings- have been reported. 
In addition to being a physical system exhibiting novel 
properties, graphene and graphene layers have been pro- 
posed as promising candidates for future nanoelectron- 
ics. The epitaxial graphene grown on SiC is of particu- 
lar interest due to the compatibility with current silicon 
technology 

Besides their unusual electronic structure^ vibra- 
tional properties and phonon spectra are also of fun- 
damental interest from which many physical properties 
(such as thermal conductivity and heat capacity) can 
be derived. Furthermore, phonons are crucial for study- 
ing the quasiparticle dynamics^ and electrical transport 
properties. Electrons excited by optical methods can 
be scattered into another states by electron-phonon cou- 
pling (EPC). It has been suggested that the scattering 
between electrons and the optical phonon modes greatly 
affects the high-field ballistic transport properties in car- 
bon nanotubesrii In graphene and metallic single-wall 
nanotubes, the EPC strongly affects the phonon frequen- 
cies, giving rise to Kohn anomalie a^^'^'^ and possible soft 
modes or Peierls distortions fi^^i^ 

Many experimental methods have been used to mea- 
sure the phonon dispersions of graphite, such as inelas- 
tic neutron scattering (INS)r^ electron-energy loss spec- 
troscopy (EELS)ri high-resolution electron energy-loss 
spectroscopy (HREELS),'IS' and inelastic x-ray scatter- 
ing (IXS).^^'^° These measurements require large enough 
samples of crystalline quality and are limited to specific 
directions or phonon modes. More recently, Mohr et al.^'^ 
have presented complete measurements of both the opti- 
cal and the acoustic phonon modes along the directions 
F-K-M-F of graphite using IXS. The results in these mea- 



surements are very close. In contrast to bulk graphite, 
Raman scattering has been widely used for probing the 
G-band in graphene layers that corresponds to the F 
phonons j^'2^'2^'^ Recently reported Raman spectra for 
graphene layers show that the intensity and position of 
the first-order G-band as well as the second order D- 
band (historically named the G' band) are modified with 
an increasing number of layers . 

On the theoretical side, Griineis et al^ presented the 
phonon dispersions of graphite using the 4th-nearest- 
neighbor force constant (4NNFC) approach. However, 
it has been argued that due to the Kohn anomaly at F 
and K, it is not possible to obtain the correct phonon dis- 
persions near F and K from the force constant method. — 
Dubay and Kressei^ performed density-functional theory 
(DFT) calculations of the phonon dispersions in graphite 
within the local density approximation (LDA). Their re- 
sults are in good agreement with phonon-measurements 
by HREELS. Using the LDA and the generalized- 
gradient approximation (GGA), Wirtz and Rubic^^ cal- 
culated the phonon dispersions of graphite and obtained 
results close to the vast majority of the experimen- 
tal data-points. At the GGA-PBE level, Mounet and 
Marzari^^ also presented a detailed calculation of the 
phonon dispersions of graphene and graphite. 

With regard to graphene layers, it is unclear how the 
phonon properties are affected by the stacking order and 
the weak interlayer coupling. This effect is important for 
understanding the EPC in multilayer graphene as well as 
the interpretation of the Raman spectra. For example, 
the phonon dispersion around K is crucial for the correct 
interpretation of the Raman second order D peak. It has 
been shown that in few-layer graphene, the electronic 
dispersions near the Fermi level exhibit various features 
depending on the stacking order. ^'^ In this work, the vi- 
brational properties of one- and few-layer graphenes are 
calculated using density-functional perturbation theory 
(DFPT).^^ The monolayer, bilayer (AB stacking), and 
trilayer (ABA and ABC stackings) are considered in or- 
der to illustrate the effects of stacking order and inter- 
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FIG. 1; Stacking structure for (a) monolayer, (b) AB bilayer, 
(c) ABA trilayer and (d) ABC trilayer graphene. 



layer coupling. The van der Waals corrections in graphite 
have been shown to be important in order to correctly de- 
scribe the long-range binding properties. ^^'"^"•■^^ However, 
previous theoretical calculations based on DFT with both 
LDAii and GGA^^ have indicated that rather reason- 
able vibrational properties of graphite can be obtained 
within DFT as compared with experiments. We find that 
the phonon dispersions for graphene and graphene layers 
exhibit somewhat different characteristics, especially at 
r and K. Detailed analysis of the phonon modes is also 
presented. 



II. COMPUTATIONAL DETAILS 

Density-functional calculations are performed us- 
ing the ESPRESSO code^^- with the LDA. TrouUier- 
Martin (TM) norm-conserving pseudopotentials^ gen- 
erated from the valence configuration of 2s^2p^ for C are 
employed. The wavefunction and the charge density are 
expanded using energy cutoffs of 110 and 440 Ryd, re- 
spectively. Methfessel-Paxton smearing^ with an energy 
width of 0.03 Ryd is adopted for the self-consistent cal- 
culations. The dynamical matrices are calculated based 
on DFPT within the linear response. For the integra- 
tion over electronic states in the calculations, we use a 
48 X 48 X 1 uniform k-point mesh. A 6 x 6 x 1 grid is 
used for the phonon calculation to obtain the dynami- 
cal matrices. We have carefully tested these parameters 
and the phonon frequencies are converged to be within 1 
cm~^. 

Figure [1] shows the two-dimensional (2D) primitive 
cells for monolayer, bilayer (AB stacking), and trilayer 
graphene (ABA and ABC stacking). A large vacuum re- 
gion of more than 10 A along the z direction is used 
to minimize the interactions between graphene layers in 
different supercells. 
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FIG. 2: (Color online) Phonon dispersions of graphite along 
the F-A direction. Solid lines are present calculational results. 
Circles are IXS data from Ref. |20, and squares are neutron 
scattering data from Ref. The dotted lines are smooth 
curves through the measured points. 



III. RESULTS AND DISCUSSIONS 



The optimized LDA lattice constant in the graphene 
plane is 2.45 A, in good agreement with the previous cal- 
culated result,^" which is also close to the experimental 
value of 2.46 A for graphite. '^^ For the bilayer and tri- 
layer systems, the lattice constant in the plane remains 
almost the same as in graphene. The optimized interlayer 
spacing is 3.33 A, slightly smaller than the experimen- 
tal value of 3.35 A in bulk graphite.— In comparison, 
we obtained a theoretical value of 3.32 A for graphite, 
which is close to previous LDA results^^S The interlayer 
binding energy (defined as the total energy difference be- 
tween the coupled and uncoupled graphene layers) of bulk 
graphite is calculated to be 25.2 meV/atom, while this 
energy falls to 12.3 meV/atom for an AB bilayer. For 
the ABA and ABC trilayer, the interlayer binding en- 
ergies are both 16.5 meV/atom. Our result of graphite 
is comparable to previous calculations using a combined 
density-functional and intermolecular perturbation the- 
ory approach.— To validate the DFPT phonon calcula- 
tions, we also calculated the low energy phonon disper- 
sions in bulk graphite along the F-A direction (perpendic- 
ular to the layers), as shown in Fig. [2] Except for a small 
frequency shift, our LDA dispersions agree with experi- 
mental data rather well. Therefore, we believe our cal- 
culations yield reliable descriptions of the phonon prop- 
erties in graphite and graphene layers. This agreement 
between LDA and experimental results indicate an error 
cancellation for energy variations near the equilibrium 
layer separation, even though the state-of-the-art local 
or semilocal exchange-correlation functionals are not able 
to properly describe the long-range interlayer interactions 
dominated by van der Waals dispersion forces (see Refs. 
[33 and m for details). 
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FIG. 4: (Color online) Pattern of atomic displacements for 
the TO/LO modes at V in monolayer graphene. 



FIG. 3: Phonon dispersions for monolayer graphene. 



A. Phonon Properties of Graphene 

Figure [3] shows the phonon dispersions for monolayer 
graphene calculated at the theoretical lattice constant, 
which will be compared with multilayer results in the 
next section. In contrast to the linear dispersion near the 
r point for the in-plane TA and LA modes, the out-of- 
plane ZA mode shows a (f' dispersion, which is a charac- 
teristic feature of the phonon dispersions in layered crys- 
tals as observed experimentally.^^-'^^'^" The same feature 
also appears in bilayer and trilayer phonon dispersions, 
as will be discussed in Section UlI Bl 

The calculated frequency (1586 cm^^) of the degen- 
erate LO and TO modes at F is slightly smaller than 
the previous value of 1595 cm~^ obtained by Dubay 
tt al.^ but is in excellent agreement with the exper- 
imental result of 1587 cm~^ by inelastic x-ray scatter- 
ing measurements.— At the Brillouin zone corner K, the 
phonon energy of the symmetric TO A[ mode (1306 
cm~^) is close to the frequency (1326 cm"-'^) calculated 
by Wirtz et ali^ Our result is also consistent with the es- 
timate by Yao et al4i from their high-voltage transport 
measurements for graphite. They suggested that scatter- 
ing by phonons with an energy of about 1300 cm^^ gives 
rise to the dramatic conductance drop at a high bias. In 
general, our calculated phonon dispersions for monolayer 
graphene are comparable with those obtained in previous 
calculationsi^i^ and agree very well with experimental 
results J^i^ 

A previous study by Piscanec et alM- showed that the 
degenerate £'29 modes at F and the highest TO mode at 
K have a strong EPC, leading to Kohn anomalies in the 
phonon dispersions. A detailed analysis of the origin of 
the strong EPC for these modes is presented below. 

For a specific phonon mode v with wave vector q, the 
displacement of atom j {j=a,(3) in unit cell m will oscil- 
late according to the following expression in the classical 
picture: 

u^J = e,Re{6^,,,,e'(^-^"--*)} 

s—x,y,z 



= ^s|e^^,Jcos(g-i?™ -cj,^t + ^^^_J, (1) 

s—x^y^z 

where Rm is the lattice vector for unit cell m, 1/)^^, ^ de- 
notes the phase factor of the complex eigenvector e^^ ^ , 
and LOqi, is the phonon frequency. 

For the degenerate TO/LO phonon modes at F, 
Figs, mja) and (b) schematically show the atomic dis- 
placements associated with the two eigenmodes. Clearly, 
two neighboring atoms vibrate opposite to one another. 
This gives rise to a large bond distortion and couples to 
electronic states near the Dirac point (which can be pro- 
jected into two states localized at atom a and /3, respec- 
tively) through an intravalley scattering (with phonon 
qwO). Therefore, a strong EPC is expected, which has 
also been demonstrated by the effective mass theory^ as 
well as the tight-binding model 42i 

In contrast, for the highest TO A[ mode at K the clas- 
sical displacements of neighboring atoms a and /3 follow 
the pattern 

u°' = uo[excos(- - ivt) - eysin(- - ujt)], (2) 
u'^ = uo[e xCos{^ - ivt) + eySm{^ - ujt)] (3) 

and atoms a and f3 move circularly. In particular, one 
moves counterclockwise, while the other clockwise, as 
shown in Fig. [5l Accordingly, each atom approaches its 
three nearest neighbors successively during one period. 
Fig. [5] shows three snapshots of the atomic displacements 
in one period. 

Since the degenerate electronic states at the Dirac 
point can be projected into two states localized at atom 
a and /3, respectively, the above mode of ionic vibra- 
tion facilitates the transition of an electron from atoms 
a to /9 or vice versa, resulting in the electronic interval- 
ley scattering via a phonon with momentum K. There- 
fore, a strong EPC is expected for this mode compared 
with other modes at Based on the detailed analyses 
of these modes, we anticipate a distinct electron-phonon 
interactions for these modes in few-layer graphene. The 
results will be presented elsewhere. 
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FIG. 5: (Color online) Three snapshots of atomic displace- 
ments for the highest TO mode at K in monolayer graphene, 
where atom a approaches its three nearest neighbors /3 (a), 
/?' (b), and /3" (c) successively. 



B. Phonon Dispersions for Graphene Layers 

In this section, we focus on the optical phonon modes 
in multilayer graphene. The phonon dispersions for bi- 
layer and trilayer are shown in Fig. [5] (a). The detailed 
dispersions for the high optical branches near F and K are 
enlarged in Figs.[6l^b) and (c), respectively. The optical 
phonon frequencies are also listed in Table HI 

Compared with the monolayer result, several distinct 
features can be identified for graphene multilayers. First, 
there is one (two) additional low-frequency mode with en- 
ergy of about 90 cm~^ at F in bilayer (trilayer) graphene. 
These modes arise from interlayer movement (so-called 
'layer breathing' modes). Second, at F the doubly de- 
generate E2g branch in the monolayer evolves into two 
(three) doubly degenerate branches for bilayer (trilayer) 
graphene, as shown in Fig. [SJb). These small splittings 
are due to the weak interlayer coupling: about 5 cm~^ 
for bilayer and no more than 5 cm~^ for trilayer (see Ta- 
ble [l|. Moving away from F, each of these degenerate 
branches breaks into two different modes. Recent exper- 
iments show that the Raman G-peak intensity enhances 
almost linearly with respect to the layer number (up 
to four layers) j ^^i^^ This phenomenon could be ascribed 
to the increased number of optical phonon modes at F 
within a small energy window for multilayer graphene. 

The stackings of graphene layers have various point 
group symmetry for the F phonons. The monolayer 
graphene possesses the De/i symmetry (Schonflies nota- 
tion). It reduces to D^d for the AB bilayer and ABC 
trilayer, and D^h for the ABA trilayer. Correspondingly, 
their high optical zone-center modes are of different mode 
symmetry: E2g mode in graphene evolves into Eg and £"„ 
for the AB bilayer, 2E'+E" for the ABA trilayer, and 
2Eg+Eu for the ABC trilayer. The Eg and E" modes 
are Raman active, £"„ is IR active, while the E' modes 
are both Raman and IR active. Therefore, a complete 
picture of the zone-center modes can be obtained from 
a combination of Raman and IR measurements. These 
mode splittings provide significant information about the 
layer number and the stacking geometry. 

In Figs. [Tlini we show the schematic atomic displace- 
ments of these optical eigenmodes at F for the AB bilayer, 
ABA trilayer, and ABC trilayer, respectively. These 



high-frequency phonons are derived from the superposi- 
tions of intralayer optical modes in each graphene plane. 
For the modes in the bilayer as shown in Fig. [3 the two 
atoms on top of each other in two adjacent layers vibrate 
either in the opposite direction {Eg mode, 1587 cm~^) 
or in the same direction (E^ mode, 1592 cm~^). Sim- 
ilar atomic displacements can also be seen in ABA and 
ABC trilayers, as shown in Figs. [5] and [D In other words, 
the original intralayer modes couple to each other via in- 
terlayer interactions, giving rise to a small splitting in 
the final frequencies. The upper and lower modes in the 
bilayer correspond to the 'in-phase' and 'out-of-phase' 
superpositions of the two intralayer modes, respectively, 
similar to the Eiu and E2g modes in bulk graphite. 

This splitting of the phonon frequencies at F can be 
illustrated using a simple model. Using the original in- 
tralayer optical modes as the basis and assuming the in- 
teraction strength between adjacent layers is e, the re- 
duced Hamiltonian for the bilayer and trilayer can be 
expressed as: 



and 



H2 = Eol + 
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respectively. Here, Eq is the energy of the intralayer 
mode, and only the first nearest-neighbor layer-layer in- 
teraction is considered. For the trilayer, a small variant 
of 6 is introduced in Eq. ^ to account for the change 
of the on-site energy in the middle layer due to the new 
geometry. (This is similar to the on-site energy variation 
due to environmental changes in electronic tight-binding 
models4^i^) Solving the secular equation det(i/— A/)=0, 
one obtains the eigenvalues and eigenvectors. 

For the bilayer, the eigenvalues are Ai^2 = EQ±e. From 
Fig. El the Eq and e can be determined: i?o=1589.5 and 
|e|=2.5 cm^^. Eq shows a small shift compared with the 
value for a single-layer graphene (1586 cm~^) as a result 
of the environmental change mentioned above. With e > 
0, the corresponding eigenvectors are 0i^2 = (liil)'^- 
This is consistent with the displacements we obtained in 
Fig. [71 The lower frequency corresponds to the out-of- 
phase superposition of the two intralayer modes (with 
respect to the motion of the two atoms on top of each 
other in two adjacent layers), while the higher one corre- 
sponds to the in-phase superposition. In both cases, the 
two intralayer modes have equal amplitudes. 

For the trilayer, the eigenvalues are Xi=Eo, X2,3=Eo + 
{S ± v^^ + 8e^)/2, with corresponding eigenvectors 
</)i-(l, 0, -1)^, 02,3=(e, {S ± VS^ + 8e2)/2, e)^. For the 
latter two eigenvectors, the ratio of the mode amplitudes 
in each layer is ai : 02 : 03=6 : {5 ± \/6'^ + 8e'^)/2 : e. 
Using the frequencies for the ABA trilayer as shown in 
Fig. [51 we obtain S Ri 3 cm~^ and e ~ 2.2 cm~^, and the 
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KM r 0.00 0.05 0.10 0.15 0.20 0.25 0.5 0.6 0.7 0.8 0.9 
r K [2ji/aJ r- k [2^^^] 



FIG. 6: Phonon dispersions for graphene multilayers. From top to bottom: AB bilayer, ABA trilayer, and ABC trilayer. 
Column (a): full phonon spectra; column (b): optical phonon dispersions near F; column (c): optical phonon dispersions near 
K. 

TABLE I: High optical phonon frequencies u) (in cm~^) at F and K for monolayer, bilayer, trilayer graphene, and bulk graphite. 
The phonon frequencies at F and K from recent DFT calculations (with LDA and GGA) as well as experimental measurements 
are listed for comparison. The point group symmetry at F (K) for the monolayer, AB bilayer, ABA trilayer, ABC trilayer, and 
graphite is Dgh [D^h], D^d (C^v), D^h {C^h), D^d {Czv), and Dg/i (£'3/1), respectively. In parentheses are the mode symmetries. 



Graphene 



AB 



ABA 



ABC 



Graphite 



Graphite exp. 



1582-', 1581" 
1588« 



K 



1586 (£29) 
1595", 1597* 
1569"^, 1581" 

1306 (A'l) 
1371", 1326" 
1289^ 1300" 
1265'' 



1587 (£9) 
1592 {Eu) 



1318 (£) 



1586 {E') 
1588 (E") 
1593 (E') 

1316 {E[, E"' 
1324 {E'2) 



1586 (Eg) 
1589 (Eu) 
1594 (Eg) 

1318 (E) 
1325 (Ai) 



1586 {E2g) 
1595 (Si„) 



1322 (E) 



"LDA, soft projector augmented wave (PAW), Ref. 

''LDA, Hard PAW, Ref. 

"LDA, TM potentials, Ref. 

"^GGA, TM potentials, Ref. 11 . 

"GGA, Ref. [il. 

•''Expt. uj{E2g), Refs.^T Egm . 

fExpt. uj{Eiu), Ref. I47l|and [50II. 

''Inelastic X-ray data of Ref. [19II and [20(1 . 



mode amplitude ratios of 2.2:5.0:2.2 and 2.2:(-2.0):2.2 for 
these two modes, respectively. This result agrees with 
the displacements from the direct first-principles calcu- 
lations shown in Fig. [5] Similar results can be obtained 
for the ABC trilayer. Based on these numerical results. 



one can easily estimate the optical phonon frequencies 
for more graphene layers using an interlayer interaction 
of 2-3 cm~^ and the values of Eq and 6 obtained above. 

The highest optical phonon branch at K becomes dou- 
bly degenerate at 1318 cm~^ in the bilayer system (see 



6 



Table m, nearly 12 cm~^ higher than that in graphene. 
The degeneracy is imposed by the symmetry of the bi- 
layer. The two degenerate modes correspond to two in- 
tralayer modes within individual layers with little cou- 
pling between them. In contrast, the optical phonons in 
the ABA split into three modes, with two of them be- 
ing almost degenerate, while the degeneracy is imposed 
by symmetry in the ABC trilayer. For the ABA and 
ABC trilayers, the highest phonon (singlet) frequencies 
are 1324 and 1325 cm~^ , respectively. This result is con- 
sistent with the Raman observation that the second-order 
D mode at about 2700 cm~^ increases with an increasing 
layer number;^ The second-order D mode in the Raman 
spectrum of graphene and graphene layers, which is dou- 
ble of the highest optical phonon frequency at K, can be 
well illustrated using a double-resonant modeli^ 

Figure [10] shows the schematic atomic displacements 
of the three eigenmodes at K for the ABA trilayer. Each 
eigenmode comprises a superposition of the intralayer Ai 
modes from each layer. As shown in Fig. (TUl the two al- 
most degenerate low-frequency modes (about 1316 cm~^) 
correspond to the combinations of the modes from the 
top and bottom A layers, while the high branch is an 
intralayer mode from the middle B layer almost exclu- 
sively. 

The splitting in frequencies at K can be analyzed in a 
similar way as before. In contrast to the F phonons, the 
coupling between adjacent intralayer modes is zero due 
to the mode symmetry, and there might only be a small 
interaction between the top and the bottom layers. In 
this case, the Hamiltonian can be expressed as: 








V 





E'o + 5' 





V 








with rj the small interaction between second nearest- 
neighbor layers. The eigenvalues are A = E'q ± rj, and 
Eq + S' . According to the frequencies as shown in Fig. [TUl 
we obtain S'=7.2 cm~^ and |77|=0.2 cm~^. The ampli- 
tudes are ai : 02 : a^—ztl : : 1 for X — Eq ± rj, and 
: 1 : for A = i?Q -I- (5'. This is consistent with the 
atomic displacements illustrated in Fig. [TO] with rj <0. 
Note that the small interactions between second nearest- 
neighbor layers have induced a small splitting of the two 
low-frequency modes. 

In the case of ABC stacking, the interlayer coupling 
matrix elements are identically zero. Therefore, the three 
eigenmodes at K are localized on each of the three lay- 
ers, respectively. Due to the higher on-site energy in the 
middle layer, one state will be higher than the other two 
modes, with the lower two modes doubly degenerate due 
to the symmetry. 

IV. SUMMARY 

In summary, we have studied the phonon dispersions 
and vibrational properties for monolayer, bilayer, and tri- 



(a) 1587 cm"' (b) 1592 cm"' 




FIG. 7: (Color online) Atomic displacements of the two split 
optical branches (a) 1587, and (b) 1592 cm~^ at F for the AB 
bilayer. Only one mode for each degenerate pair is shown. 



(a) 1 586 cm'' (b) 1 588 cm"' (c) 1 593 cm 




FIG. 8: (Color online) Atomic displacements of the three split 
optical branches (a) 1586, (b) 1588, and (c) 1593 cm~^ at F 
for the ABA trilayer. The length of the arrow represents 
the amplitude of the eigenvector. Only one mode for each 
degenerate pair is shown. 



layer graphene using the density-functional perturbation 
theory. Due to the weak coupling between layers, the 
highest optical phonon branch at F in graphene splits 
into two (three) doubly degenerate branches with small 
yet unnegligible splittings for bilayer (trilayer) graphene. 
The splitting is about 5 cm^^ for the bilayer. In trilayer 
ABA and ABC graphene, these splittings are about 2 
cm~^ and 5 cm~^, respectively, which are not equally 
spaced. These modes are of various mode symmetry and 
exhibit different sensitivity to either Raman or IR mea- 
surements and therefore a combination of Raman and IR 



(a) 1 586 cm ' (b) 1 589 cm ' (c) 1 594 cm 




FIG. 9: (Color online) Atomic displacements of the three split 
optical branches (a) 1586, (b) 1589, and (c) 1594 cm"^ at F 
for the ABC trilayer. The length of the arrow represents 
the amplitude of the eigenvector. Only one mode for each 
degenerate pair is shown. 
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(a) 1316.2 cm"' (b) 1316.5 cm"' (c) 1324cm"' 




FIG. 10: (Color online) Atomic displacements of the three 
optical phonon modes (a) 1316.2, (b) 1316.5, and (c) 1324 
cm"'^ at K for the ABA trilayer. The length of the arrow 
represents the amplitude of the eigenvector. Only one mode 
for each degenerate pair is shown. 

measurements of the zone-center optical modes should 
give a clear identification of the layer number as well as 
the stacking geometry. 

A simple interaction model is applied to illustrate the 
frequency splitting and the characteristics of the eigen- 
modes at T. The interlayer coupling strength is identified 
as about 2 cm"-'^. In the trilayer system, a shift of about 
3 cm^^ in the on-site energy in the middle layer is deter- 
mined. 



The frequency of the highest optical phonon mode at 
K in bilayer (trilayer) graphene is about 12 (18) cm~^ 
higher than that in monolayer graphene. For trilayer 
graphene, the K-A'j^ mode splits into three branches in 
the ABA trilayer, with the two lower modes nearly dou- 
bly degenerate. It is found that the on-site energy varia- 
tions for the middle layer in ABA and ABC are about 7-8 
cm~^, higher than that of T phonons. Due to the sym- 
metry, the interlayer coupling between adjacent layers for 
these intralayer modes is zero. 
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